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Abstract 

We consider the M/G/\ queue with a processor sharing server. We 
study the conditional sojourn time distribution, conditioned on the 
customer's service requirement, as well as the unconditional distribu- 
tion, in various asymptotic limits. These include large time and/or 
large service request, and heavy traffic, where the arrival rate is only 
slightly less than the service rate. Our results demonstrate the possible 
tail behaviors of the unconditional distribution, which was previously 
known in the cases G = M and G = D (where it is purely exponen- 
tial). We assume that the service density decays at least exponentially 
fast. We use various methods for the asymptotic expansion of inte- 
grals, such as the Laplace and saddle point methods. 



*Department of Mathematics, Statistics, and Computer Science, University of Illi- 
nois at Chicago, 851 South Morgan (M/C 249), Chicago, IL 60607-7045, USA. Email: 
qzhen2@uic.edu. 

'''Department of Mathematics, Statistics, and Computer Science, University of Illi- 
nois at Chicago, 851 South Morgan (M/C 249), Chicago, IL 60607-7045, USA. Email: 
knessl@uic.edu. 

Acknowledgement: This work was partly supported by NSF grant DMS 05-03745. 



1 



1 Introduction 



One of the most interesting service disciphnes in queueing theory is that of 
processor sharing (PS). Here every customer in the system gets an equal 
fraction of the server or processor, and this has the advantage that shorter 
jobs get served in less time than, say, under the first-in-first-out (FIFO) 
discipline. 

PS queues were introduced in the 1960's by Kleinrock (see [I], [2]) and have 
been the subject of much research over the past 40^ years. In these models 
one of the main measures of performance is a given (also called tagged) 
customer's sojourn time distribution, conditioned on that customer's service 
time. The sojourn time is the time the tagged customer leaves the system 
after being served, assuming the customer arrives at time zero. 

We denote by V(x) the conditional sojourn time, with x being the service 
time. If the tagged customer arrived to an empty system and no further 
arrivals occurred in the time interval [0,x], then V(a;) = x. But in general 
V(x) > X as the tagged customer must share the server. We denote by b{x) 
the service time density, by p{t\x) the conditional sojourn time density, and 
by p{t) = b{x) p{t\x)dx the unconditional sojourn time density. We note 
that p{t\x) has in general a probability mass along t = x, but p{t) is generally 
continuous. 

The M/M/l-PS queue assumes Poisson arrivals and i.i.d. service times 
with density b{x) = fie~^^. In [3], Coffman, Muntz and Trotter obtained an 
expression for the Laplace transform oi p{t\x) (i.e., for £[6"*^^"^^], where s is 
the Laplace transform variable). In [3] Morrison removed the conditioning on 
X and studied p{t) in the heavy traffic limit, where the Poisson arrival rate A 
is nearly equal to the service rate fj, (thus p = A//i t 1). Setting e = 1 — p, in 
[1] asymptotic results were obtained for the time scales t = 0(1), t = 0(e~^) 
and t = 0{e-^). Most the mass is concentrated in the range t = 0(e ), 
and the asymptotic series involves modified Bessel functions. For moderate 
traffic intensities with p < 1, the tail behavior of p{t) is given by 

p(t) ~Cor^/^e-^«*e-^«*'^', t^oo (1.1) 

where = p (1 — ^/p)'^ and the constants Bq and Cq are given in |5]-[7]. The 
result in (11. II) was obtained for the M/M/1 queue under a random order of 
service (ROS) discipline, but there is a close connection between the waiting 
time distribution in the ROS model and the sojourn time distribution in the 
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PS model. This relation, along with some extensions, is explored in [S] and 
[H]. In [TU] we studied the conditional density p(t\x) for various asymptotic 
ranges of x and t, for both a fixed p < 1 and in the heavy traffic limit where 
Ptl. 

A more difficult model is the M/G/l-PS queue, where the service density 
is general. This was analyzed by Yashkov in [11], [12] and by Ott [13j. These 
authors obtained an explicit, albeit complicated, expression for £[6"*^*^^-*]. 
Inverting the Laplace transform leads to an expression for p{t\x) as a contour 
integral (see fl2.4p ). but the integrand is a nonlinear function of another 
contour integral, which is in turn defined in terms of the Laplace transform 
of the service density. In the case of deterministic service times, where b{x) = 

— l//i), much more explicit results are available (see [13] and [S]). Also, 
the tail behavior of the unconditional sojourn time density was derived by 
Egorova, Zwart and Boxma p3j as 

p(t)~C"e-^'* {G = D), t^oo (1.2) 

where the constants A' and C are explicitly characterized in [H]. Comparing 
f fO) and ([L2l) we see that the tail behaviors of the M/M/l-PS and M/D/l- 
PS models are quite different. 

In this paper we will study both the conditional sojourn time density 
p{t\x) and the unconditional density p(t) in the M/G/l-PS model. As in [10] 
we shall consider various asymptotic limits, such as x and/or t — >■ cxo with a 
fixed p < I, and 1 — p = e — * 0+ with space and time scaled by e. Here p = 
Ami where mi = xb{x)dx is the mean service time. We consider service 
densities b{x) that have "thin tails" , with decay that is at least exponential as 
X — > oo. We shall show that the basic asymptotic structure of the conditional 
density p(t\x) is essentially independent of the service density (though the 
formulas do depend on the Laplace transform of b{x)). In contrast, the 
unconditional density is highly dependent on the tail behavior of b{x). We 
shall make specific assumptions on this tail, first assuming that 

b{x) ~ Mx'?e-^^\ r > 1, (1.3) 

where M, (> 0) and q are constants. Thus (11. 3p allows for roughly ex- 
ponential or even thinner tails, such as a Gaussian. Then we shall discuss 
"zero-tail" service densities, by assuming that b{x) has support for < a; < A 
(e.g., b{x) = 1/A corresponds to uniformly distributed service times). In the 
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zero tail case the behavior of p{t) as t — oo and p < 1 is determined by the 
behavior of b{x) near the upper hmit of its support, and we will assume that 

b{x) ^ a, {A - xy'-\ xTA (1-4) 

where a,, and i/* are positive constants. 

We will obtain a wide variety of tail behaviors of p{t) as t ^ oo for the 
general M/G/l-PS model, that are different from either (11 .11) or (11.21) . We 
shall also identify the class of service densities that lead to purely exponential 
tails, such as G = D in (II. 2p . 

We mention some related work on various PS models. Ramaswami [15] 
studied the G/M/l-PS queue and obtained explicit results for the uncondi- 
tional moments of the sojourn time. Various asymptotic properties of the 
conditional and unconditional moments and distribution were derived in [16] . 
The G/G/l-PS model has not been analyzed exactly, but some approxima- 
tions are discussed in Sengupta [17] and the tail exponent of the unconditional 
sojourn time density was derived by Mandjes and Zwart [TS]. Specifically, in 
[TB] the authors characterized the limit \og[p(t)] —Aq ast —>■ oo, assum- 
ing that the arrival and service densities have at least exponential tails. In 
[19j Zwart and Boxma analyze the M/G/l-PS queue with heavy tails, where 
the service density has algebraic or sub-exponential behavior as x ^ oo (thus 
AT = 0, r = 0, g < -1 in ([L3D, or < r < 1). 

For PS models one is also interested in the sojourn time conditioned 
on the number of other customers in the system when the tagged customer 
arrives. The conditional sojourn time for the M/M/l-PS model, conditioned 
on this number rather than the service time x, was studied by Sengupta and 
Jagerman [20] and Guillemin and Boyer plj. A good recent survey of sojourn 
time asymptotics in PS queues is in Borst, Niifiez-Queija and Zwart [22] . 

In this paper the main methods used are for the asymptotic expansion of 
integrals, such as the Laplace and saddle point methods, and good general 
references are the books of Bleistein and Handelsman [53] and of Wong [21] . 

The remainder of the paper is organized as follows. In Section 2 we 
summarize and briefly discuss our main results (see Theorems 2.1-2.5). In 
Section 3 we derive the results for p{t\x), for moderate traffic intensities 
p < 1. In Section 4 we consider p(t\x) for p I 1, and various scalings of space 
and time. We remove the condition on x in Section 5, treating both p < 1 
and p ~ 1, and here we make the assumptions (II. 3p or (11.41) . 
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2 Summary of results 



We assume that customers arrive according to a Poisson process with rate A, 
at a single processor-sharing server. The customers' random service requests 
are i.i.d. random variables with density function b{y), and Laplace-Stieltjes 
transform 6(r) = e^^^ b{y)dy. We assume that 6(r) is an analytic function 
of T for 3?(r) > — cq for some cq > 0. Thus all the moments of the service 
time are finite, and we set 

POO 

mu= yH{y)dy={-lfh^^\T)\r=, {ken). (2.1) 

^0 

In particular we denote the service rate by ix where 

- = mi= yb{y)dy = -b'{0). 
^^ Jo 

Then the traffic intensity is defined as p = A//i and we assume that p < 1 . 

The Laplace transform of the equilibrium sojourn time distribution, con- 
ditioned on the tagged customer (or job) requiring x units of service, was 
derived by Ott \13l and Yashkov [11], [12], who obtained 

prp-sv(x)i ^ ^ - P (2 n) 

^ J (l-p)G^is,x) + sGsis,xy 



where 



' -rx^ , ^ , r-A(l-6(r)) 
e Gi{s,x)dx = 



and 



r[r-s-A(l-6(r))] 

re-G.is,x)dx^ P[^-M1-K^))] 
Jo r2 [r - s - A (1 - 6(r))] 

Thus, the Laplace transform of the denominator in the right-hand side of 
(ED is 



/(r;s) = / [il-p)G,is,x) + sGsis,x)]e-^^dx 

[1 - p) - (1 - p) A (1 - 6(r)) r + spr-sA(l- 6(r)) 



r2 
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(2.3) 



Taking the inverse Laplace transform of fl2.3l) . fl2.2l) becomes 



E[e-«v(x)| 



l-p 



(2.4) 



where Br^- is a vertical contour in the complex r-plane, on which 3?(r) > 0. 

By taking the inverse Laplace transform of (12.41) . the probability density 
of the sojourn time, conditioned on service time x, is 



p{t\x) 



2m J ^ J 



l-p 

2m 



Br, 



2m 



— e^-f{r-s)dT 



Err 



ds. 



(2.5) 



Here Brg is a vertical contour in the complex s-plane with 9ft(s) > 0. Note 
that p{t\x) will in general have a probability mass along t = x. Analyzing the 
integral (12. 5p . we obtain the following expansions for p{t\x), valid on different 
space and time scales. 

Theorem 2.1 For a fixed p < 1, the conditional sojourn time density has 
the following asymptotic expansions: 

1. X ^ oo, t — X ^ with x{t — xY = 0(1), assuming that 
b{y) a y"'^ , as y — > (a, z/>0), 



p{t\x) 



1-p 
2m JB,. 



XaT(u) 



X 



ds 



{l-p)5{t-x) e-^^ 

— ^ m,' V um.) 



um—1 



m=l 



2. x,t — > oo, 1 < t/x < oo, 

, {l-p)rl[l + \y{n)]"' 
p(t\x) ~ 



Sl \ 27TxXb"(T' 



(2.6) 



(2.7) 



where = s*(t/x) and = T^{t/x) satisfy the equations 

t - X 

-b\n) = J e-^*yyb{y)dy=^, = - X {1 - b{r,)) . (2.8) 
3. x,t^ oo, t/x^ = 0(1), 



P(t\^) 93/2 ^1/2 . 2 ,5/2 ^ " ^ ^ >^ exp ' 



/2(2n+l)2x2 
X ^ 2t 



n=0 

X^ \ 

2tl, 

A6"(^o) 

where Sq = s*(oo) and tq = r*(oo) satisfy: 

X 6'(ro) = -1, So = To - A (1 - 6(ro)). (2.9) 

4- X = 0(1), t ^ oo, 

p(t|a;) ~ (1-p) J(x)e^=(^)*, (2.10) 



/ f{T;s)e^^dr) 

JBrr J 

and Sc[x)(< 0) is the maximal real solution of 



J{x) = ^U- 

ds V 2712 



(2.11) 



S = Sc{x) 



/ /(r;s,)e-^dr = (2.12) 



or 

l — p + xSc + —^. I — c/r = 0. 



2m 



Br. r2 [r-s,-X{l-h{T))\ 



The result in Case 4 was also recently derived by Yashkov [25], who 
characterized Sc{x) and J(a;) in a different form. 

In the asymptotic matching region between Cases 3 and 4, we have 



p{t\x) ~ (1 p)7r^AVoMyXTo))' jsot + + Ot/xn 

2sqX'^ 



(2.13) 
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where 



2 ' 3ro • I • J 

Then result (12.131) holds for x,t ^ oo with t = 0{x^). It can be extended to 
larger ranges of t, e.g., to t = 0{x^), by including an additional factor of the 
form exp{Dt/x'^). 

We note that in Case 2, iit/x ~ l/(l-p), by (ES]) we have s,(l/(l-p)) = 
0, r,(l/(l-p)) =0 and 



2 Xm2X 



X 



1-p 



+ 0({t 



X 



Then the formula fl2.7p simplifies to the Gaussian 



p{t\x) 



P) 



3/2 



V2 TT A m2 X 



exp 



2 A m2 X 



1-p 



1-p 



which gives the spread about the well known mean value E [V(x)] = a;/(l— p). 

As a special case of the M/G/l-PS model, we consider the M/i?fc/l-PS 
model, in which the Erlang service time density function is given by 



h{y)= ^ 1 (A:GN,y>0) 



and thus 



6(r) 



k jj, 



k fi + T ^ 

Then we obtain the following more explicit results. 



(2.15) 
(2.16) 



Corollary 2.1 For the M/Ek/l-PS model with traffic intensity p < 1, the 
conditional sojourn time density has the following asymptotic expansions: 



1. X oo, t — X ^ with x{t — xY = 0(1), 



p{t\x) 



1-P 
2m 



Brs 



XikfxY 



X 



ds 



k—l „~\x 



(l-p)5(t-x)e-^^' + ii^^(^^) ^('-^^'^ ' 



(k-iy. 

X oF,([];[l + i,l + p...,2-l,2];Ap'=x(t-x)'=). (2.17) 
Here oF/^d ]; [bi, 62, bk]', z) is the generalized hypergeometric function. 
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2. x,t ^ oo, 1 < t/x < oo, 



p{t\x) 



where 



I 1 fc+2 

(1 - p) kfl [t - {pkY+^ it - x)]Jii (A/i'^)— (^) (f ) 



tsl v/27r(A: + l) 



X 



= (A /)'=+! 



1 - x/t 



k + l 



(2.18) 



(A: + l-x/t)-A-/iA;, (2.19) 



f^: = kp 



P \ k + l 



_Vl - x/t 



(2.20) 



5. x,t^ oo, t/x^ = 0(1), 



p{t\x) 



{1- p)k^k{k + l){l- p^^)^ 
2y/27FJI p^oh) [p- (A; + l)piTT +A;]^t5/2 

(2n + l)2A;(A/i^)^x2 



n=0 
2fc 



2(A; + l)t 



X 



^ X - 2 1 



where 



so = (A; + l)(A/)^-(A;/i + A), 
fo = k {X p'')'kTi - k p. 



4- X = 0(1), t — > oo, 



p(t|x) 



(1 -p) e"'-(^)* 



(2.21) 

(2.22) 
(2.23) 

(2.24) 



where Ti = rj(s) (z = 1, ...,k + 1) are i/ie A; + 1 poles of f{r] s) in l{2.3\) . 
with residues Ri{s) = Res{f,T = ri{s)), and Sc{x) is the maximal root 
o/E?=ie-«W-i?,(s) = 0. 
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In the matching region between Cases 3 and 4, we have 

(1-p) {k+lfTx\l- p^^Y 



p{t\x) 



2 pfe+i [p- {k + 1) + kY 

X exp j/c// (1 - pfcTT)x + sot + Bt/x'^ + Ct/x^^ , 



(2.25) 



where 



TT^ (A; + 1) 



2 /c/zpfc+i 
(A; + 1) [{k - 4) p^i - {k + 2)] 

S/c^yU^pfeTT (pfeTT _ l) 



(2.26) 
(2.27) 



We note that for the M/Ek/l-PS model, u = k, a = {k p)^/{k - 1)!, and 
r*, s*, So and tq in Theorem 2.1 are explicitly computable, as given by f*, s*. 
So and tq in Corollary 2.1. 

We next consider the heavy traffic case, where A t A*- Letting e = 1 — p 
(thus e — > 0"*"), we have the following results for general service time distri- 
butions. 

Theorem 2.2 For p = 1 — e, where e 0"*", we let t = T/e and x = X/e. 
The conditional sojourn time density of M/G/l-PS model has the following 
asymptotic expansions: 



1. x^O{l),t^O{l), 



St 



2. x = 0(1), T = 0(1), 
p{t\x) 

X (x 



e^^[T- p + ph{T)] 
27ri JBrr t^[t- p + p b{T) - s] 



dr 



ds. 

(2.28) 



~ — e ' — e 



5(T) (T - 2x) e-^/^ 



X^ 



2-ni 



x^ 



C_ T^[T-p{l-b{T))] 



dr 



(2.29) 
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where 
Q*{x) 



90 + 90 1712 ^3 + 15 m2 (4 mg — 3 m2 ^4) x 



270 /im^ 

+ (9 + 20 rrig — 30 m2 ms 7714) 



(2.30) 



and TUi is the i*'* moment of the service time distribution, given by 
l[2.1\) . The contour C_ can be taken as the imaginary axis in the t- 
plane, indented to the left of t = Q, where the integrand has a pole of 
order 4- 

3. X,T = 0(1), T - X 0+ with T-X = T^ e^+^Z" = Oie^+^Z"), assum- 
ing that 

b{y) ~ a y'^^^, as y {a,iy > 0), 



p{t\x) 



^-nX/e^l~l/u 

27ri 



a J. r/ia;r(z/) ■ 



dS 



Brs 

m=l 



.(2.31) 



4. X = 0{1),T = 0(1), and 1 < T/X < 00, 

T s^- X 



p{t\x) 



,3/2^2 



exp 



+ T(f* - 



'2 71 fiTb"{%) 
where = s^{T / X) and = f^(T / X) satisfy 

1 + fib'{f^) = X/T, = - |U (1 - b{f^)). 



(2.32) 



(2.33) 



5. X = ^Z = 0{y^),T = 0{l), 
p[t\x) ~ 

V/im2 vr T 



^exp 



n=0 



(2^+1)2^2 

2 yum2 T 



(2.34) 



6. X = 0(1), T = 0/e = 0(e "'^), we g'Zi'e the expansion in three different 
forms: 
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(a) 



p{t\x] 



e' / X 

~ : exp 



e 



B., (l + V?)^-(l-v^)^exp(-|;|ye) 



(2.35) 



(b) 



p(t\x) ^Vexp 

(2n)! 



2(n+ 1)X _ 4 
4 



2?i 

xE(-i) 



Z!(2n-0! 



(2.36) 



X 



/im2 



i:'_i_i(2;„) +i:'_i_2(^r, 



ifere Di,(-) is the parabolic cylinder function and = ''^^"^^.^e ■ 



(c) 



p{t\x) e^^e'''^''-^'^G{vn). 



(2.37) 



n=l 



ifere Vn = Vn{X) are the real positive roots of the equation 



exp 



2X 



IV 



1 + \ 2 
1 — iv 



(2.38) 



G{v„ 



2^nexp(^) 



(vl - 1)X cos (^) + [2vnX + {vl + 1) /.m^] sin (^) ' 

(2.39) 



s<i(^;„) = Bi{vn) + Ci(^;„) e + O(e'), 



where 



2 /xm2 ' 



(2.40) 



(2.41) 
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Cl{Vn) 



K + 1: 



+ (3 /im^ t;^ - 3 mg t;^ - 3 /im^ + mg) xj . (2.42) 



For very large times, corresponding to ^ 1 (thus t ^ e ^), we have 



p{t\x) 



2 „sdivi)e 



e e 



(2.43) 



where Vi = Vi{X) is the unique root of f l2.38p in the interval (0, /im2 vr/X). 
Then the first term in the sum in (12.371) dominates. 

For the M/ E^/l-PS model, we again get more explicit expressions. 

Corollary 2.2 For the M/Ek/l-PS model in heavy traffic, we have the fol- 
lowing expansions of the conditional sojourn time density. 



1. x = 0(1), t = 0(1), 



zm J Br. ^ 



2. x = 0(1), T = 0(1), 



e™[(r-/i)(A;^ + r)'= + /i(A;^)^] 



-dr 



-1 



ds. 



(2.44) 



a;2 5{T) + (T - 2x) e"^/^ 



X'' 



k+1 

X ( g*(x) + ^ 

i=3 



Qj + kfi Q ^ 



where 



Q*{x) 



{k + 1) Q| 



(2.45) 



270(.+ l).Vn^°^^^)^^^ + 9°^^ + 2H^^)^^^ 
+ 15A; /i (A;2 + A; - 2) X - (A;^ + 9k^ + 6k - 16) 



and Qj , j = 3, k + 1 are nonzero roots of 

{Q - ^i){Q + k ^if + ^i{k ^if = 

Note that Q = is a double root and the other roots have ^{Q) < 0. 
For example, if k = 2, = —3 fi; if k = 3, Q3 = (—4 + i^/2) /i and 
Q4 = (-4 - zV2) /i. 
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3. X,T^ T-X ^0+ withT-X^n e^+^/^ = 0{e^+^/^), 



p{t\x) 



Q rp 

e * exp 



Brs 



— s'^ — J 

{k-l)\ 



X 



oF,([ ]; [1 + i, 1 + I 2 - i 2]; Z+^XTj^) . 



^. X = 0(1), T = 0(1) and 1 < T/X < oo, 



p{t\x) 







f ^ ^ 2{fe+i) 


4fc+3 

v/27r(A; + 1)X (| j ^(^+^' 


;(^)^(^ + l)-(fc + l)]' 



X exp < — ijl 



X exp /X 
5. X^y/^Z^O{y/-e),T^O{l), 
p{t\x) ~ e'/2 



8k 



{k + l)TTT 



E 

n=0 



exp 



{2n+iykfiZ^ 
2{k + l)T 



6. X — 0(1), T = 9/e = 0(e we have the following three different 
forms of the expansion: 



(a) 



p{t\x) 



e'^ Ilk 



{k + l)'Ki 



- exp 



fikX nkQ 



X 



k + 1 2{k + l) 

( iJ.k0 
I 2(fc+l) 



exp f - f^Vt) 



(1 + VC)^ - (1 - VI? exp ( - ^x/e) 
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(b) 



p{t\x) 



X 



n=0 
2'+3/2 / 



exp 



k + 1 



4 



i+3 



1=0 



(2n)! 



fee 



/! {2n-iy. 



where z„ 



k ^ 



(fc+i)e 



3[e + (2n + l)X]. 



n=l 



i/ere f„ = fn(X) are the real positive roots of the equation 



exp 



2kfiX 

—, —iv 

k + l 



1 I 2 
1 + 



1 — iv 



(2.46) 



where 



CAVn) 



2VnX + 



;^-l)Xcos(^) + ^ 



(fc+l) (vl+l) 



sm 



/ k ^iv„ X \ 
\ k+l ) 



BAVr, 



kfl{dl + l) 



{2k + 1) 



2{k + l) ' 

2{k-l)kiiXvl 
~ k^i{vl + l)X + 2{k + l) 



6(fc + l)2 

For very large times with ^ 1, we again have 

p{t\x) ~ e^e'-^^^'^^G^Vi), 

where Vi = Vi{X) is the unique root of fl2.46p in the interval (0, ■^^^-^)- 

The asymptotics of the conditional sojourn time distribution are generally 
dependent on the service density only through equations such as (12. 8p or 
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(12.91) . which involve the Laplace transform. This is true for all the scales in 
Theorem 2.1 except Case 1, where t ^ x and the expansion depends on the 
local behavior of the service density b{y) as y 0~^. We now examine the 
unconditional sojourn time density, defined by 

p{t) = [ h{x)p(t\x)dx. (2.47) 
Jo 

The structure of p{t) is highly dependent on the behavior of h{x) for x — >■ 
oo, and we need to make specific assumptions on the tail of the service 
distribution. We first assume that the service time density function behaves 
as 

h{y) ~ My''e-^y\ y ^ oo {M, N > 0, 1 < r < 2). (2.48) 

We remove the condition on x by using the results in (I2.13P (if p < 1) 
or (12.431) (if p ~ 1), and evaluate asymptotically the integral in (I2.47p . We 
thus derive the following results for the unconditional sojourn time density 
as t — s> oo, from which we realize the variety of possible tail behaviors for the 
M/G/l-PS model. 



Theorem 2.3 As t ^ oo, the unconditional sojourn time density has the 
following asymptotic expansions, assuming the tail behavior in ^2.4^ , and 
that p is fixed with p < 1. 

1. Ifr = l, 

p(t) ~ait^exp(sot-7it^/^), (2.49) 



where 



M(l -p)r2 42+5 ^ 22+1 

V6s^{N + ro) 3 

r (iV + ro)[6 6"(ro) + ror(ro)] 
X exp < 

I 3ro6"(ro) 
and Sq and Tq are computed from ^2.9\) . 
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2. // 1 < r < 2, 



pit) ~ a2 exp(so t - (f){^,,t)). 



Here 



= {Nc- ^) + roe^^ - ^t^: 



(2.50) 
(2.51) 



(1 — p)Afrn ,1^ 3 . 2(9-1) y,, , ,-, 

a2 = _^ Q -(A^r)r + 2 71-2+ r + 2 [A6 (Tq)] 



3 I q-1 
'2~'~ r + 2 



where B and C are given by \2.1Ji\) , and = ^*(t) satisfies 



IB \ 3C 
0^(e,t) = [NrC~^ + —)t^^ +Tot^^ +—t'^^ =0. (2.52) 



VKe have the following expansion(s) of and (f){C^.,t) as t oo 



Nr 



r+2 i_ 



{r + 2){Nr 

.2 1^2 \ til 



4 

I r+2 



tr+2 



+ < 



(5-r).„M^'A6"Ml^^^_ ^^^^ 

2(r + 2)2(Arr)^ 
6 6"(ro) + ror(ro) 



[ (r + 2)ro6"(ro) 



t '■+2, 



3/2 < r < 2, 



0(e*,t) 



!l±^(Arr)42 [vr^AS'Yro) 
2 r L 



+ To 



Nr 



r + 2 



'^^(^^0-^pA6"(ro) 
7r2A[6 6"(ro)+ro6'"(ro)] 



2-r 

r+2 2-r 
tr + 2 



tr+2 

1 < r < 3/2, 



3 To 



r+2 

t^, 3/2<r<2. 



5. //r = 2, 

where 



L7r2A6"(ro)J 

p(t) ~ as exp(so t - 73 1^/' - ^3 1'^'), (2.53) 

(3ArC-roS)2 ND 



2+5 
4 



(l-p)M7r''/2+ir2[A6"(ro)] 
«3 = .a+5 „ exp 



16NB^ 



B 
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73 = 7ry^2iVA6"(ro), 
NC + tqB 

^~ iVl/4|5|3/4' 



and 



D 




For the M/ E^/l-PS model, we note that the parameters in the tail of the 
service time density in (12.481) are M = {kfi)''/{k — 1)!, N = k^, q = k — 1 
and r = 1. Thus Case 1 of Theorem 2.3 applies, with ai and 71 given by 



and So is given by sq in (I2.22p . 

We next consider the unconditional distribution in the heavy traffic limit. 
We find that the expansion of p{t) in (I2.47P is different for the time scales 
t = 0(1), t = O(e-i) and t = 0(e"^-2). For t = 0(1) we can use fICTD 
to evaluate the integral in ( 12.47p . but no further simplification is possible. 
Below we give the results on the large times scales, and we note that most 
of the probability mass occurs where t = T/e = 0(e~^). 

Theorem 2.4 For p = 1 — e, where e 0^ , the unconditional sojourn time 
density has the following asymptotic expansions. 



ai = 




11 = 1^'^' [fc(fc + l)]V3(/A)m), 



1. t = T/e = O(e-i) 




(2.54) 



18 



where S{T) is given by the double integral 

1 



S{T) = - I h{x) 







2tii 



Uc^ r2 [r-/.(l-6(r))] 



-dr 



dx. 



Here the contour C+ is taken as the imaginary axis in the r-plane, 
indented to the right o/r = 0. 

2. t = (T*/e"+2 = 0(e""-2) With l<r<2, 

V2 7r /ima M X'^ G {v^ (X) ) e^/^"? 



Pit) 



a, {v[{X))^ + Vi{X) <(X)J + r (r - 1) ^im^N X-'' 
X exp [-7/;o(X(a,), a,) e"" + Ci(t;i(X)) a, e^"^] . (2.55) 

i/ere X = X{a^:) satisfies 

_ N r X^-^ [/im2 X {vj (X) + 1) + 2 /x^ mj] 

.?(x)K(x) + i] ' ^'-''^ 

ipo{X,a^:) is defined by 

MX,a,) = N X' - B,{v,{X)) a,, 

Vi = Vi{X) is the unique root of Ii2. 38\) in the interval (0, /im2 vr/X), 
and Bi{vi) and Ci(fi) are given by i2.41\) and (2^J^- 



3. t = a/e^ = 0(e"4) with r = 2, 

V2 TT /ima M G{vi (X) ) e^"'? 



[K(X))2 + t;i(X)<(X)] +2/im2X 
X exp I - MXia), a) t'^ + Ci(t;i(X)) a + D^{vi{X)) a 

4/im2Ar + 2cr[K(X))2 + i;i(X)<(X)] /■ 
ifere X = X(cr) satisfies 

2NX [/im2 X (i;2(X) + 1) + 2 /i^ m^] 



vl{X) K(X) + 1] 
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and Di{v) is defined by 



DAv) 



72 [X (t;2 + 1) + 2 /ims] ^ 



where 



dQ = 2AjjL — 12 jj, m2 + 8/im2 m3 + (f ^ + 1) m2 m4 — (w +3)7713 



Aj/ ml 



36 ii' {v' - 3) - 24 /i {Av^ - 3) 777,2 "^3 



+3 (3 t;^ - 2 77^ + 3) 7772 7774 - (11 -A2v^ + 27) 777^ 



2 (11712 



36 / (f ^ + 3) 777^ - 24 /i (t;^ + 2 1;^ - 3) 777^ 7773 



+ 3 (3t;^ - 7v^ - Tv"^ + 3)77737774 - (13 - 27t;^ -A5v'^ + 27) 777^ 



(t> + 1) — 36 /i 7772 + 24 /i7772 7773 — 3 (t; + 1) 777 2 7774 



-(5t;^-18 7;2 + 9)777^ 



If the service density had even thinner tails, say with r > 2 in fl2.48p . 
we can easily extend Theorem 2.3. The main complication is that we would 
need further terms in the expansion of Sc{x) as x ^ 00, which has the form 
So + B/x'^ + C/x^ + D/x'^ + 0{x~^). When r = 2 the 0{x~^) term affects 
the leading term in the expansion of p{t) as t —>■ 00. The asymptotic eval- 
uation of (12.471) involves balancing the factors exp(i?t/a;^) and exp(— A^x^), 
which occurs when x = 0(^'^^). We can also extend Theorem 2.3 to more 
complicated tail behaviors of b{y) of the form 

b{y) ^ M y" - N y'- + Niy'^ + ■ ■ ■ + Ni y''] 

where < < r;„i < ■ ■ ■ < ri < r. This would be needed, for example, if 
b{y) were a truncated Gaussian centered at some non-zero y (then r = 2 and 
ri = 1). We shall next consider densities with "zero-tail", and these lead to 
different behaviors of p{t). 
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Now we assume that the service time density function b{y) has finite 
support foT < y < A and behaves as 

biy)^a,{A-yy'-\ y ] A, («„ z/, > 0). (2.58) 

The structure of the conditional sojourn time density p{t\x) is the same as in 
Theorem 2.1, but the unconditional sojourn time density p{t) is determined 
by the behavior of b{y) near the upper limit A of its support. 

As t — > oo with fixed p and p < 1, we remove the condition on x by using 
the results in fl2.10p . For the heavy traffic case with p ~ 1, we remove the 
condition on x by using fl2.29p on the large time scale t = 0(e~^). We thus 
have the following results for the unconditional sojourn time density. 

Theorem 2.5 The unconditional sojourn time density has the following asymp- 
totic expansions, assuming the service density behavior in Ii2.58\) . 

1. t oo with p fixed and p < 1: 

where J{x) is given by ^2.11\) and Sc{x) is the maximal real solution of 
(EJl). 

2. t = T/t = 0{t-^), p=l-e with e ^ 0+.- 

j-A -Tlx 

p(t) = e b(x) dx + 0(e'^). (2.60) 

Jo X 

For very large times with T ^ oo, fl2.60p becomes 

p{t) ^ ta,V{u,) A^'^'-^T--' e-'^/''. (2.61) 

We can show that when fl2.59p is expanded in the heavy traffic limit p | 1 , we 
also obtain fl2.6ip . We note that as p j 1, (12.121) shows that Sc{A) ~ —e/A = 
0(e) and thus s'^{A)t ~ T/A^. The 0{e^) term in (IZeOD is the same as that 
in f l2.54p . except that the integral over x is truncated aX x = A. 

If b{A) is non-zero and finite then z/* = 1 and f l2.59p shows that p{t) has 
an exponential tail with the additional algebraic factor 1/t. This additional 
factor disappears only in the limit of z/* — > 0, but then (12.580 shows that b{y) 
develops a probability mass aX y = A. Our results show that p{t) will have 
a purely exponential tail only if b{y) consists of one (then G = D) or several 
point masses, or if b{y) has a point mass at the maximum of its support 
{y = A) with all the remaining mass in the range < y < A^, with A^ < A. 
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3 Brief derivations of the conditional sojourn 
time density for the case p < 1 



In this section, we first give a brief derivation of the conditional sojourn time 
density for the M/Ek/l-PS model, with a fixed traffic intensity p less than 
one. Then we sketch the derivations for the general service time density. 

The Erlang service time density and its Laplace transform are given by 
f l2.15p and fl2.16p . Then (12.31) has the following explicit form: 

f{r;s) 

_ r(l - p) [{t - A) (r + fc /i)^ + A(fc ^u)^] + sp [(r - /i) (r + fc /x)^ + fijk /i)^] 

~ r2 [{r-s- A)(r + k p)'' + X{k p)>'] ' 

(3.1) 

From (13. ip we see that the numerator of /(r; s) has a double zero at r = 0, 
so / is analytic at r = and the poles r = r(s) satisfy 

{t - s - X){t + k p)'' + X{k p)'' = 0. (3.2) 
Consider first the limit x,t oo with 1 < t/x < oo. We define F by 

F(s,x) = -^/ /(r;s)e^^cir. (3.3) 

From (13.21) we have 

By applying the residue theorem in (13.31) and using (13.41) . we have 

F(s,x) ~ i?(r(s),s)e"(^)^ (3.5) 

where 

p/ / N N (r(g) + kp) 

" T^{s) [{k + 1) r(.) -ks + k{p~X)] ^"^-^^ 

and r = r(s) is the largest real root of (13. 2p . The other poles of / lead to 
exponentially smaller terms and (13. 5p holds when x —>■ oo and 3fJ(r) > on 
the vertical contour i?r^. We next define 

ip{s) = ip(s;-) = s- -t{s), (3.7) 

\ X/ X 
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and then by fl2.5l) we have 

The integrand in (13.81) has saddle points where ^p'{s) = so that there is 
a saddle point along the real axis at s*, where 5* and = r(s*) are the 
solutions of 

(r - s - A) (r + kfi)'' + \{k^if = 0. 
From (13.21) we have 

{k + l)r{s) + k{f,-s-Xy ^^-^^^ 

Solving the system ^M) with the help of flXTOj) leads to fl219|l and ^Ml- If 
we shift in (13. 8p to Br'^, on which 3fJ(s) = and use the saddle point 
method (see, e.g., Wong [24]) with the steepest decent direction arg(s — s*) = 
±7r/2, we get 

p{t\x) ~ , e"^(^'*). (3.11) 

A/2 7ra;(/9"(s^,) i?(r(s^,), s^,) 

But, from ( 13. 7p and ( 13. 2p we have 

/// N ///X A: (fc + 1) (r(5) + kfi){s- t{s) + A) 
09 Is) = —T Is) = o — . (3.12) 

[{k + l)T{s) + k{fi-S-\)Y 

Using (EJ]), (EZD and fl312D in (13111) leads to (I2A8D . 

Next we consider x, t ^ cxd but with t/x ~ 1. The previous calculation 
is not valid since, from (12.19p . the saddle point s^: oo. From (13.21) . f^, has 
the following expansion as s,,, — cxd: 

r, = s, + A ^f-^ + O ' 



■'jl. ^ s* 



We return to (13.81) and note that i?(r(s^,), s,,,) ~ 1 as s* — ^ cxd. Then we ap- 
proximate the integrand for s large (more precisely we can scale s = 0{x^^'') 
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with x{t—xY = 0(1)) to get the integral representation in fl2.17p . Expanding 
exp[\{k fi)''s~~''x] as a geometric series, using the identity 



1 

27ri 



6(t — x) m = 0, 
(km — 1)! ~ ' 



and the generahzed hypergeometric function 



E 



'^.F.(ll;|H-il + f,....2 4,21;.), 



'^{m+l)\{km + k-l)\ {k - 1) 

we obtain the second expression in fl2.17p . 

Now we consider x, t — oo but with x/t small. From (12.191) we let x/t — >■ 
0, then the saddle point 5* sq, which is given by (I2.22p . Then from (I2.20p . 
has the following expansion 

= To ± fa a/ S-Sq + 0{S - So), 

where tq is given by (12.231) and 



Ta 



2k{Xii^)—^ 



k + l 



This means that on this scale, we must re-examine (13.30 as now two poles 
of the function /(r; s) determine the asymptotics of the integral (13. 3p . We 
denote these two poles as f_|_ = f_|_(so) and f_ = f_(so)- Then for s ^ Sq 



F{s,x) ~ R{f+,So)e^+^ + R{f.,So)e^- 

Rq 



Vs - So 



(3.13) 



where R{r±,so 
Rq = ■ 



=^ and Rq is given by 



1 



k^k + 1) 



k + l 



p - {k + 1) pfc+i + k 



1-p 



k + l 



(3.14) 
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We return to fl2.5l) . shift the contour Br^ to Br", which is shghtly to the 
right of So, and use fl3.13p . thus obtaining 



p{t\x) 



2m .1 B^ 



i t — TQ X 



V-s - So 



// Ro (e 



TaVs-SQ X 



faVs-SQ 



ds 



1-p 

Rn 



-TO X 



OO ^ 

y— 

^ 2m 



n=0 



X / V -5 - So exp 



st — (2n + 1) Ta Y s — So X ds. (3.15) 



Using the identity 



^ . , V'5 - So exp 

1 r(2n + l)2f2a;2 ^ 2 i 



s t — (2n + 1) Tq a/s — So X 



c/s 



exp 



(2n + l)2 f2a;2. 
At 



Sat- 



in flXT^ . we obtain ([MI]). 

Finally, we consider the 0(1) and t ^ OO. Now all the A; + 1 

poles of / in fl3.ll) contribute to the asymptotics. We denote these poles 
by Xj = rj(s), i = 1, ■ ■ ■ , k + 1, which are the solutions of (13. 2p . and (13.31) 
evaluates to 



fc+i 



F(s,a;) = ^e^"W^i?,(s), 



(3.16) 



i=l 



where Ri{s) = i?(rj(s),s) is the residue of /(t; s) at r = rj(s). Then (12.51) 
becomes 



p{t\x) 



1-p 



1 



-ds. 



(3.17) 



J Br. E':=ie^^^'^^Ri{s) 

From (I3.17P we obtain (12.240 by locating the pole s = Sc{x) with the largest 
real part, which is the maximal solution of F{s,x) = 0. 

We can simplify (I2.24p for x — oo, which leads to an explicit expression 
in the matching region between Cases 3 and 4. As x — oo, we have Sc sq 
and the expansion 

~ ~ ^ B C ^(\ 

Sc = So + - + ^ + ^ + — 
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We then expand fl3.2p about r = tq, to find that A = 0, B < (which 
corresponds to Sc < Sq) and that r has the following expansion: 



To ± 21 



2kjip^ ~ 1 
k+1 ' ' X 



+ 



k + 2 
_3{k + l) 



k ^ p'' 



k + l 



2{k + l) \B\ 



i + ofi). (3.18) 



Then as in the analysis of Case 3, two poles of the function /(r; s) dominate 
the expansion of F(s,x). We denote these two conjugate poles as ri and 
T2 = Ti. From (13.61) we have 



R^{s) = R{n{s),s) 



HP" 



k + l 



2k^{k + l)\B\ 



p — {k + 1) pfc+i + k 



1-p 



k+l 



X 



{k-A)pT^ -{k + 2) p-{k + l)pT^ + k 
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pp" 



k + l 



3k^k + l) {1-p—^y 
p-{k + l)p^ + A;"'^ 



{2k\B\)n 



1-p 



k + l 



+ 0{- 

^x 



(3.19) 



and R2{s) = Ri{s). Thus (I3.16P is approximately 

F{s,x) ~ e"^W"i?i(s) + e"^(^)"i?2(s) 
~ 2K(e"i(^)^i?i(s)). 



(3.20) 



Expanding (13.201) as x — > oo and noting that Sc is a root of F{s, x) = in 
this limit, we obtain the expressions for B and C in (12.261) and (12.271) . Then 
we can simplify (I2.24p to 



p{t\x) 



(1 - p) exp {sot + B t/x'^ + C t/x^) 



d_ 

ds 



But, 



2 3fJ(e-i(^)^i?i(s)) 

R!^{s) + Ri{s)xr[{s) 
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(3.21) 



(3.22) 



and 

R[{S,)=Oix^), (3.23) 



1 



r,(..)~-.^^^^x. (3.24) 



Using (Km . (1321 and (Km . (Km becomes 



as V 



,2 



2 

k + 1 



p - {k + 1) pfc+i + k 



1 

1 — pfe+i 



2 



which leads to (Km . 

We now consider general service densities b{y). The basic scales in The- 
orem 2.1 are the same as those for the case in Corollary 2.1, but some 
of the definiting equations are more complicated, becoming transcendental 
rather than algebraic. We consider the function (12.31) and note that 6(0) = 1 
and b'{0) = — l//i, so r = is not a pole of /(r; s). The poles r = r(s) of 
/(r; s) now satisfy 

r-s- A(l-6(r)) =0. (3.25) 

For the case x, t — > oo, with 1 <t/x < oo, the asymptotics are obtained 
analogously to the Ek case and (13. 7p still applies, but (13.61) . (I3.10p and (I3.12p 
must be replaced by 

(3.26) 

and 

A6"(r) 



(1 + A6'(r))'' 

The first equation in (12. 8p follows from using ip'{s) = 0, (13. 7p and (I3.27p . 

For the case x,t —>■ oo with t/a; ~ 1, we have to make some assumptions 
about the behavior of b{y) as y ^ 0. We assume that b{y) ~ ay'^~^ (a, u > 
0) for y ^ 0^ . Then the Laplace transform of the service time satisfies 
6(r) ~ ar{i')T^'^ as r ^ oo. Then (12. 6p is obtained in the same way as 
(I2.17p . although we cannot express the sum in (12. 6p as a hypergeometric 
function for non-integer u. 
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Next we consider x,t oo but with x/t small. From fl2.8l) . letting 
t/x — * oo, we have Sq = s*(oo) and tq = r^,(oo), which are given by fl2.9p . 
Then t^, has the following expansion as s — > sq: 



A6"(ro) 



+ 0{s-so) 



Thus again two poles of /(r; s) dominate the expansion of F{s,x) and the 
calculation is similar to the Erlang case, with fl3.14p becoming 



Rn 



2A&"(ro) 



Finally, we consider the = 0(1) and t — > oo. For the general service 

time distribution, all the singularities of the function /(r; s) contribute to 
F{s,x). Then fl2.10p is obtained by using the residue theorem at the largest 
pole Sc{x) of the integrand in (12. 5p . which is the maximal real solution of 
F{s,x) = 0. 

In the asymptotic matching region between Cases 3 and 4, we let x ^ cxo 
and 

ABC 

Sc = So + - + — + — + 0[ — ), 

X X x^ \x / 
and expand (13.250 at r = tq. We find that A = 0, B < and two conjugate 
poles Ti and T2 of the function /(r; s) dominate the behavior of F{s,x). 
Analogously to (I3.18P and (13.190 . ti and Ri{s) have the following expansions: 



Ti = ro + 



2|5| 1 r \B\b"'{To) 
A 6" (to) X ^ [3A(6"(ro))2 



iC 



2\b"{To) |SH 



- + 0(4), 



Ri(s) 



+ 



I Sr 



X — 



TiJ2\V'{r,)\B\ 



I sic 



2BTi 



2A6"(ro) \B\ 



si [&h"{T,) + T,y"{To 
3Ar3(6"(ro))' 
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and T2 = Ti, R2{s) = Ri{s). The constants B and C are obtained by expand- 
ing fl3.20p as a; — s> oo using F{sc{x),x) = 0, and this leads to fl2.14p . Using 
(K2T^ with (§0, B, C) replaced by (sq, B, C) and 



X 



d 



ds V 

we obtain (12.131) 



vrA6"(ro) 
So 



vrAro6"(ro) 



4 Brief derivations of the conditional sojourn 
time density for the case p ~ 1 

Now we consider the M/G/l-PS model with a traffic intensity that is close 
to one, and let p = 1 — e with < e <^ 1. 

First, we consider x = 0(1) and t = 0(1). Using A = /x + 0(e) we obtain 
from ([23D 

/(^; s) = — + 0(e). 

r2 [r-s-/i(l-6(r))] 

This leads to (12.281) . On this scale the solution does not simplify much, but 
there is little probability mass in heavy traffic on the time scale t = 0(1). 

Next, we consider x = 0(1) but for large time scales t = T/e = 0(e^^). 
In (13.31) we replace p as 1 — e and scale s a.s ew, and we have 



F{s,x) = ^ / e 



T X 

Brr 



T + W w'^ 2 ^/ S\ 

T- e + — ^ T e + ^ e 

r2[r-p(l-Kr))] 



dr 



e{l + wx) + e^w^^ [ — -dT + 0{e') 

27rWB,^ r2 [r-/i(l-6(r))] 
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Then from (12. 5p . we obtain 



p{t\x) 



-,wT 



27rz Jbt^ F{ew,x] 



dw 



27ii 
e 

—e 

X 



-,wT 



ew 



Br^ [l+WX 2m{l + Wx)^ Jsr, T^[t - fl{l-b{T))] 



-dr 



T/x _ ^2 



6{T) ^ (r-2x)e-^/" 



x^ 



x^ 



27ii jBrr r2 [r - - 6(r))] 



dw 

dr. 
(4.1) 



The function 



g{T,x) 



r2 [r-^(l-b{T))] 
has a pole at r = of order 4. By the residue theorem we have 

77— I 9{^,x)dT = ReSr=o{g{T,x)) + ^ I g{T,x)dT. 
2vr« JBrr 27rz Jc_ 



(4.2) 



Here we shifted the contour Bvj. to C_, which can be taken as the imaginary 
axis in the r-plane, indented to the left of r = 0. Then we define 



1 d^ 

Q,(x) = ReSr=o{g{T,x)) = — lim — 



^2 gTX 



r-/i(l-6(r))J 



(4.3) 



which leads to fl2.30p . Note that we assumed that all the moments of the 
service time are finite, which are given by (12. ip . Expression (12.290 is obtained 
by using (14. ip . (14. 2 p and (14. 3p . The term proportional to 6(T) in (14. ip does 
not mean that there is actually mass at T = 0, but rather corresponds to 
the small (0(e)) mass that exists in the shorter time scale t, where (12.280 
applies. 

Now consider x = X/e = 0{e-^) and t = T/e = 0{e-^) with 1 < T/X < 
oo. By the same argument as in Section 3, the pole r = r(s) of /(r; s) with 
the largest real part satisfies (I3.25p . We replace A by /x (1 — e) in (I3.25p . which 
yields 

r-s-fxil-e){l-biT))=0, (4.4) 

and then expand t as r = Ta + r^e + O(e^). Then = Ta{s) and = Tb{s) 
satisfy 

Ta - S - fl {1 - b{Ta)) =0 
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and 



n 



/i(Kra)-l) 



1 + ;u6'(r„) 

In (13.71) we replace t and x by T/e and X/e respectively, to get 



T 
T 

s — -Tais)-nis)e + Oie^) 
Ms) + fi{s) e + 0{e^) 



and rewrite (I3.8p as 



p{t\x) ~ — 



exp [Xipo{s)/e + X(pi{s 
R{r{s),s) 



ds. 



(4.5) 



(4.6) 



Here R is as in (I3.26p . with A replaced by /x. Then the integrand in (14. 6 p has 
a saddle point where v^o(s) = 0, which satisfies 



T T 1 

X X l + ^6'(r,) 



Ta - S - yU (1 - b{Ta)) = 0. 



0, 



(4.7) 



We denote the solution of (14.70 (T/X) and = s^{T/X)^ which 

leads to (I2.33p . Then by the standard saddle point method, (14. 6 p asymptot- 
ically evaluates to 



p{t\x) 



,3/2 



iy/2 vrX v9q(s*) R{%,s^[ 
But by (13:261) and (03]), we have 



exp 



Xipo{s*)/(^ + Xipi{s^) 



(4.8) 



and 



¥?o(s,) 



f2[l + /x6'(f, 

/ife"(f. 



[1 + /i6'('r* 



3 • 



(4.9) 



(4.10) 
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Using dH]), fHTTU]) and (USD in (gSD, we obtain ([232]) • We note that if 
X ^ oo and T ^ oo but T/X = 0(1), the approximation fl2.32l) remains 
vahd. 

For the case X = 0(1), T = 0(1) and T - X 
that the service time density behaves as b{y) ^ ay 
We note that the saddle point s,, oo as T/X - 
find that r has the following expansion for s — * oo 



0+, we again assume 
''-^ (a, z/ > 0) for y ^ 0. 
1. Then from (14.41) . we 



s + /i 



/xa r(z/) 



Following the same argument as in Section 3, we can easily obtain fl2.3ip . 
once we scale s as 5e^/^ and let = {T - X) e-^-^l"" = 0(1). 

Next, we consider X = = O(v^), T = 0(1). If we let T/X oo 

in fl2.33p . it follows that f^, — > and the saddle point — > 0. By (14.41) and 
scaling s = 0(e) = ew, we find that r has the following expansion: 



± 



2w I- wm^ 
Ve + 



3 7712 



/X7772 



3 



e. 



Now two poles of the function /(r; s) dominate the behavior of F{s,x). We 
approximate F{s, x) by the sum of the residue at these two poles, where from 
(I3.26P we also have i?(r(s), s) ~ ±a/ fim2we/8. Then from (12.50 we obtain 



p{t\x) 



exp 



riwZ 



Br-w ^Jw 1^1 — exp ( — 

rp OO 

^exp 



W 



n=0 



2V2wZ 
(277+ 1 



dw 



f^wZ 



dw, 



where the contour Br^ is a vertical line in the ty-plane slightly to the right 
of 777 = 0. Then (12.341) follows by using the identity 



Br^ 



exp 



(271 + 1)^2^^ 



dw 



exp 



(277 + 1)2^2 

2 fj,m2 T 



We note that by using the Poisson summation formula 



oo 

E 



oo oo 



^2iriym 



i^{y) dy, 



32 



where \l/ is the Fourier transform of ip, we can rewrite fl2.34p as 



p{t\x) 



,3/2 



oo 

l + 25^(-irexp( 



n=l 



n 71 fim2 T 
2Z2 



(4.11) 



From (14. lip , we can easily verify that Cases 2 and 5 in Theorem 2.2 asymptot- 
ically match, in the intermediate limit where x — ^ oo and Z 0. Similarly, 
Cases 4 and 5 match in the limit where X — and Z ^ oo, which follows 
easily from (12.341) . 

Now we consider X = 0(1) and T = 6/e = 0(e^^) (thus x = 0(e^^) and 
t = 0(e"^)). Similarly to the previous time scale, two poles, at ri = ti{s) 
and T2 = T2(s), dominate the behavior of F{s,x) and we have 



p{t\x) 



27xi J ReSr=TAfir; s) e^^) + ReSr=T2ifiT; s) e^' 

ds. 



ds 



27ri J Br, Ri{s) e^i(^)^ + i?2(s) e^^C^); 



(4.12) 



We scale s = O(e^) by setting s = (y^;^) and then from (14. 4p and (13.261) 
we have 



n,2 



/im2 



e and i?i,2(s) = R{ti,2) ~ ±^i^-^^-e. 



which leads to (12.351) . 

Furthermore, we expand the integrand in (12.351) as a geometric series, and 
we have 



p{t\x) 



\ 



/im2 7n 



( ^ 
exp 

V /im2 2 yum2 



exp 



0< 



2 fj.m2 



Br, (1 + V^)^ 



n=0 ^ 



2n 



exp 



{2n + 1)X 
/im2 



(4.13) 



Note that if we let X ^ oo, then the n = term in (14.130 dominates, and 
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we have 
1 



exp 



271"« J Br. (1 + V^)^ L2/im2 /^"^2 



7r0 



+ 



26 2^7772 + X + e 



erfc 



2/i7772 



7r0 

Here we used 



exp 



7ryU7772 



2/777i20 



/im2 

X ^ oo. 



V^/2/ir7720 



exp 



X2 



2/77712 6 



2 r ,2 , 1 



erfc(z) = 1 — erf (2;) 



Then (14.131) becomes, for fixed and X — > 00, 



23/2^2 



: exp 



X 



e 



e , as z 00. 



X2 



(4.14) 



^/¥jIrrl^ 1/77772 2/i7r72 2/ir7i26- 

When X = ©(e^^) and T = 0{e-^) but with T/X = 0(1), (l2:32|) remains 
vahd, and letting T/X ^ 00 in (I2.32p regains (14.141) . This again verifies that 
these two cases asymptotically match. 

We return to (14.131) . and let a/^ = z — 1, with which the integral becomes 



2^ exp 

n=0 



e ^(277 +l)X 



X 



2 111112 /77r72 
(2-1)2(^-2)2" 



^2n+2 



exp 



e 2 _ + (2^ + 1)^ ■ 

-2/i7r72 /t7r72 



dz. 



Here the contour C+ can be taken as the imaginary axis in the 2;-plane, 
indented to the right of 2; = 0. Using the binomial expansion 

.2n _ (-1)^2^)! 



f l4T3l) leads to 



z-2) 



00 2n 



j=0 



{2n) 



j\{2n-j)\ 



p(t\x) W , . 

/^"^2£^^j!(2r7-j)! 



2)^ exp 



2(77+ 1)X 



X 



27r7 



exp 



/ e 



(4.15) 



-2/777i2 



^2 — AnZ )dz, 
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where 



and 



A. 



e + (2n + l)X 
2/im2 



1 



K = 2n- i + 2. 

We express the integral in fl4.15p in terms of parabohc cyhnder functions, 
using 

thus obtaining 

p(t\x) - 



Br 



27r^m2 ^ 



^exp 



2in + l)X zl 
fim2 4 



2n 

i=o 



(-lV(2n)! 
j!(2n-j)! 



;-2) 



K 



i 



e \ 



/im2 



Aim2 



(4.16) 



where 



Q + {2n + 1)X 



Replacing 2n - j by /, KTB leads to ^Ml- 

If we let X — > and 6^0 with X/a/B (thus Zn) fixed, the term with 
/ = in fl2.37p dominates and we have 



p{t\x) 



23/2^2 °° 



23/2^2 



^exp 



n=0 



E 

n=0 



exp 



2(n+ 1)X _ zl 

111712 4 - 

(2n + l)^X^ 
2;um20 



^o(2;„) 



(4.17) 



Here we note that X^/B = Z'^/T and used the fact that Dq{w) = e"'^^/'^. 
Since (14.171) is the same as fl2.34p . we have shown that Case 5 is really a 
special case of Case 6 in Theorem 2.2. 

Alternately, we can treat the problem on the (X, 0) scale by evaluating 
(14. 12p using the residues of the integrand at all the poles Sp, which satisfy 



Resr=r,is,){fir; s) e"^) + i?es,=,,(,^)(/(r; s) e"^) = 0. 



(4.18) 
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We let Sp = Aie + Bie^ + de^ + O(e^). Then from and flX^ we find 
that Ai = and we have the following expansions: 



/im2 



where 



/im2 

2m3 — 3fim2 + fim2m3Bi 

S/i^mgCi + 3/im2 — 2m3 + 3fim2{fiml — m^)Bi 
S/i^mgVl + 2/im2i?i ' 

i?i,2(s) ~ (i±r„)e±rfee^ 

_ 1 + /im2-Bi 
~ 2V1 + 2/im25T' 

[2/im2Ci — (Sfiml — m3)i?i] 
6(1 + 2/im25i)3/2 ■ 
By expanding the left-hand side of (14. 181) about e = 0, we obtain 

2r, + l 



and 



where 



n 



-2tcX 



2r„-l 



and 



r, = 2(r2-l/4)reX 



(4.19) 



(4.20) 

Setting a/1 + 2/im2-Bi = m + w (f > 0) we find that all the roots of (14.191) 
are on the imaginary axis, and with -u = (14.191) becomes (I2.38p . Denoting 
the m}^ positive solution by f„ = t>„(X), we obtain Bi = Bi{vn) by (I2.4ip . 
By solving (I4.20p for Ci, we have 



3fim2{fim2 — m^)XB'l + [(3/im| — 2m3)X + yum2(3/im| — m^)\Bi 



which leads to (I2.42p with the help of (12.411) . Note that f„ = Vn{X) have the 
following asymptotic expansions: 



/im2 2{iim2Y ^ 4(/im2)^ 



X 



X' 



X3 



n-n + 0[^^, X ^ oo, 



(4.21) 
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and 



{n — l)7Tfim2 2 4 



X ' {n — 1)71 {n — l)^7r^yum2 

Now fl4.12p becomes 



X + 0{X^), n>2,X^0. 



But 



and 



p{t\x) ~ 



1 ^[e^^^'^'^^ils) + e^2Wxi?2(s) 



/ s=s„ L 



,ri(sp)a:: 



1 .vl-1 



- + 



4t;„ 



^iispj 1 



2^ 



(1 + <)/im2 
4 f 3 e 



Thus we obtain 
d 



ds 



^ e 4^2 



S = Sp(Vn) 

2vnX + nm2{vl + 1) 



4^2 



exp 



— 1 + iv 
/im2 



X 



and 



^fe^^(^)^i?2(s; 
ds V 



s=Sp(t;„) 

Using fl4:22|) and g^l, we define G by 



as V 



s=Sp(i;„) 



23? 



(4.22) 



(4.23) 



4f2 



4^2 



exp 



— 1 + iv 
jjm2 



X 
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which leads to fl2.39l) . and then we obtain fl2.37l) with Sd(f„) = Sp(f„)/e^. 

If we consider even larger time scales, with B ^ 1 (thus t ^ e^^), then 
the largest pole Sp{vi) dominates. Here Vi = Vi{X) is the unique root in the 
interval (0,/im27r/X) of (12:381) . This leads to (^M>- The expression (12:371) 
with (I2.40p applies for time scales up to 9 = 0(e~^) (t = 0{e~^)), but for 
even larger time scales we may need further term in (I2.40p . e.g., the O(e^) 
correction to Sp. We will discuss this more in Section 5. 

The M/ Ek/l-PS results in Corollary 2.2 follows from Theorem 2.2 by 
using the j^^ moment 

(fc + j-1)! 

5 Brief derivations of the unconditional so- 
journ time density 

The structure of the unconditional sojourn time density is highly dependent 
on the tail behavior of the service density. First we assume the service time 
density function behaves as (II. 3p or ( I2.48p . For p fixed and less than one, the 
major contribution to the integral in (I2.47P will come from the asymptotic 
matching region between the scales x = 0(1) and x = 0{\/i), with t — >• oo. 
In this region, the conditional sojourn time density is given asymptotically 
by 

For 1 < r < 2, using (12.130 and (I2.48P in (12.470 . the unconditional sojourn 
time density behaves asymptotically as 

p(t)r^aoe'°* x'^~''^exp( -tqx-Nx'' + ^ + ^)dx, (5.1) 
Jq V x^ x'^ / 

{l-p)M'K^\^Ti[h"{T,)Y 



where 



2 si 



Scaling x = ^t^+^ = 0(1^+^), (15.10 becomes 



r+l 
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where t) is given by (12.51 1) . Here 5 > so as to avoid integration through 
^ = in the case of g < 3, and 5^1. By using the Laplace method with the 
major contribution coming from C,* = C,*{t), which satisfies 0^ = or fl2.52p . 
(15.21) becomes 

Mt)~^^e-^t^e--^«-'*), (5.3) 



where 

r fiR 



r 12C r-1 
fr + 2 



QB 



r{r -1)NC - —^t—. (5.4) 

If r = 1, then by (^M), ^* = {^Y^^- Using (^Ml and with 
and r = 1 in ([53]), we obtain i^^. 

If 1 < r < 2, then by (12.521) the leading term in the asymptotic 

expansion of satisfies Nr^^~^ + |f = 0, which leads to = (^l^')'^^- 
Then we can rewrite (15.31) as 

Pit) ~ eg-^t^ 

which leads to (I2.50p . We give three terms of asymptotic expansion for 
and 0(iC*, t) in Theorem 2.3, as t —>■ oo. We note that the third terms in these 
expansions are different according asl<r<3/2orr>3/2. 

If r = 2, the above analysis is still valid but we need to include the 
additional factor exp(Dt/x^) in (I2.13p . and then in (15. ip . The constant 
D is obtained by refining the approximation (12.131) so that it applies for 
X = 0{t^/^). Thus (E2]) and ([53D become 

p{t) ~ aot'^e'^' e"'exp(- + 0(e,t))c?e 

~ /, ,N ^o ^ ^ exp fsot-0(^.,t) + -4 ). 
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Here = {\B\/Nf/^ and 

A = ^0 + ^ t-4 + o(r 2), 

Thus, after simplification, we obtain (12.531) . 

Now we consider the unconditional distribution in the heavy traffic limit, 
again assuming that the service time density function behaves as (12.481) . 

For the time scale t = T/e = 0(e~^), we use (I2.29p in (12.471) . which leads 
to (12.541) after we integrate from x = to x = oo. 

To compute the unconditional density p{t) on the scale t = cr^,/e^'^^ = 
0(e-'-2) with 1 < r < 2, we use fICT) and fl2^ in (HJ^ with 6 = a./e'". 
Scaling x = X/e = 0(e^^), (12.471) becomes, since et — ^ oo, 

P(t) ~ ^ r G{v,)X'^exp\ - ^^{X,aMx, (5.5) 
Jo '- ^ 

where 

^(X,(T,) = NX"- - Bi{vi{X))cx, - Ci{vi{X))a,e 
= i!o{X,a^) +'4>i{X,a^)e. 

Hence (15.51) is a Laplace type integral, and the major contribution will come 
from where ip is minimal, which should satisfy 

^MX,a,) = NrX^-' - A_B,{v,{X))a, = 0. (5.6) 

But from (ICTj) and (l238|) . we have 



j^B,{v,{X)) - ^^^x[vl{X) + l] + 2i^^ml- ^^'^^ 



Using (15.71) in (15. 6p . we obtain (12.561) . This defines X = X(o"*) implicitly. 
Denoting the right-hand side of (I2.56P as ^1{X), we can verify that Q'{X) > 0, 
so that Q{X) is a monotonically increasing function. As we discussed in 
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Section 4, Vi{X) ~ yum2vr/X as X — > oo and vi(X) ~ ^J^l^mi^^JX as X — 0"*". 
Then fi(X) ^ oo as X ^ oo and fi(X) ^ as X ^ 0+. Hence there is a 
unique positive root of r2(X) = cr^,, which we denote by X = X(cr^,). Then 
we use the standard Laplace method in flS.Sp to get, for 1 < r < 2, 



~ , ^ = exp [- V'(X,a,)e J, (5.8) 

Y^V'o(^,t^*) 

which leads to (12.551) . 

If r = 2, to compute the unconditional sojourn time density on the time 
scale t = cr/e^ = 0(e~^), we need to include the Di{vi)e^ = O(e^) term in 
Srf(fi) in (I2.43p . -Di(fi) is obtained in the same way that -Bi(fi) and Ci{vi) 
are derived, which we discussed in Section 4. Analogously to (15. 5p and (15. 8p . 
with a^: replaced by a and r = 2, we have 

Pit) 



Jo 


oo 

G{vi)X''e 




V 


'2^M T 
e 


X'^ G{vi{X)) 


/ 92 

yaT- 


MX, a) 


X exp 1 




-' + D^{v^iX))a + 



(5.9) 

where X = X(ct) satisfies (I2.56p with a* replaced by a and r = 2. This leads 
to (1^371) . 

Next we assume that the service time density function b{y) has finite 
support for < y < A and behaves as (I2.58P near the maximum of its 
support. As if: ^ oo with fixed p < 1, we remove the condition on x by using 
the results in (I2.10p . The main contribution comes from x = A, and we have 

.A 

p{t) ~ / a,{A-xy'-\l-p)J{x)e'^^^^'dx 
Jo 

~ (1-p) a, J{A) e^=(^)* / [A - xY'~^e''^^^~^ (5.10) 



Setting X = A- u/{s'J^A) t), flHTTOl) becomes 



, , (l-p)«, J(A)e«=(^)* /-^ 
p(t) ~ ^ ^-TTT-TT-r^ / e "c?M. 
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Using u''*-^e-''du = r(z/*), we obtain fl239D . 

For the heavy traffic case, we remove the condition on x by using fl2.29p . 
and on the large time scale t = 0{e~^) we obtain (12.601) . For even larger times 
with T — >■ oo, by using the Laplace method (with the main contribution from 
X = A), (12.601) becomes 

pit) ~ '-^^ / {A - xr-' e'^^-^^l^'dx. 

This leads to ^IMh . 



A Appendix 

We will give a brief derivation of the Laplace transform of the conditional 
sojourn time distribution with deterministic service density b{y) = 6{y — 
This was derived by Ott (see (5.16) in [13]) and more recently in 
[H]. However, these authors use arguments that are specific to the case 
G = D. Here we point out that these results also follow easily from the 
general M/G/l-PS model. 

We rewrite (5.16) in [13] as 

E[e-v(V,.)| = {l-p){X + sfe-'-"^ 

' ' s^ + \[s + (\-p)(\ + s)]e-i^'li'' ^ ' 

where we replaced z in [T3] by 1. To prove flA.ip . in view of that (12. 4p we 
need to prove that 



2ni Jsr^ (A + s^^ 



But, by (O, 



/(r;.)^^i^4^ + 



r2 r^[r-s-X{l-e--/^)] 
and we have 



27!-i ./r,_ r2 



^'^ = 1 - p + s/jj.. 

BTt 
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Thus (1A.2P is equivalent to proving the following identity: 

1 r e^/^ 



27ri JBrr r2 [r - s - A (1 



e 



■dr 



_^ (l-p)A2 + (2-p)sA s + /i(l-p) 



(A.3) 



(A + S)2 (A + S)2s2 

If we scale t = fiT and set w = p + flA.31) becomes 



dT=--^. (A.4) 



We shift the contour Btt to the right so that 3fJ(T) > 3fJ(w). Then upon 
expanding the integrand in flA.4p as a geometric series and multiplying flA.4p 
by p^ we must show that 



oo 



2« T2 (T - - ^2 • (^■^) 

For L > 1 we can close the integration contour in the right half of the T- 
plane, and these integrals all evaluate to zero. For L = we close in the 
left half-plane, where there is a simple pole at T = w and a double pole at 
T = 0. Calculating the residues leads to ( 1A.5I) . thus proving ( lA.ll) . 
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